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ABSTRACT: Inspired by Mott's (1929) analysis of particle tracks in a cloud chamber, 
we consider a simple model for quantum cosmology which includes, in the total Hamil- 
tonian, model detectors registering whether or not the system, at any stage in its en- 
tire history, passes through a series of regions in configuration space. We thus derive a 
variety of well-defined formulas for the probabilities for trajectories associated with the 
solutions to the Wheeler-DeWitt equation. The probability distribution is peaked about 
classical trajectories in configuration space. The "measured" wave functions still satisfy 
the Wheeler-DeWitt equation, except for small corrections due to the disturbance of the 
measuring device. With modified boundary conditions, the measurement amplitudes es- 
sentially agree with an earlier result of Hartle derived on rather different grounds. In the 
special case where the system is a collection of harmonic oscillators, the interpretation of 
the results is aided by the introduction of "timeless" coherent states - eigenstates of the 
Hamiltonian which are concentrated about entire classical trajectories. 



1. INTRODUCTION 



The focus of attention in quantum cosmology is the Wheeler-DeWitt equation, 



= 



(1.1) 



Here, the wave function is a functional of the gravitational and matter fields on a three- 
surface, and it describes the quantum state of a closed cosmological model [1]. The most 
striking and conceptually problematic aspect of this equation is that it does not involve 
time explicitly, severely complicating efforts to extract predictions from it [2,3]. Amongst 
the many attempts to understand this feature, one is to claim that "time", and indeed 
entire histories of the universe, are already contained amongst the arguments of the wave 
function, hence no time label is required [4,5]. Whilst these claims seem to be true at 
some level in simple models of quantum cosmology, it presents us with the interesting 
challenge of reformulating standard quantum theory without the explicit use of time, and 
then demonstrating the emergence of time and of classical trajectories. Although the 
Wheeler-DeWitt equation in the form (1.1) is unlikely to be the last word in quantum 
gravity, it does seem likely that whatever replaces it will still be of a timeless nature. The 
loop variables programme of Ashtekar and others, for example, certainly preserves this 
feature [6]. It is therefore of interest to investigate this feature in simple models. 

Many attempts to use and make sense of Eq.(l.l) have been made. These attempts 
focus on simple (minisuperspace) models, in which one has an A'"-dimensional configuration 
space, C with coordinates x, and the Hamiltonian operator has the form 



The signature of the metric is typically hyperbolic so the Wheeler-DeWitt equation is like a 
Klein-Gordon equation in curved space with a spacetime dependent mass term. Associated 




-V2 + F(x) 



(1.2) 
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with it is a Klein-Gordon current 



J = z(*V** -**V*) (1.3) 

Like the Klein-Gordon equation, however, this does not produce a positive probability 
density except in very special cases (namely when there is a Killing vector associated with 
H). It also vanishes for real wave functions. One might also consider the Schrodinger inner 
product, 

(*l|*2) = j/^'' //(x)*|(x)*2(x) (1.4) 

where /x(x) is an appropriate measure, but the norm (\E'|^') typically diverges. In practice, 
most uses of the Wheeler-DeWitt equation rely on something like the "WKB interpreta- 
tion" , in which in the oscillatory regime the wave function is written in the form ^ = Ce^'^, 
where 5' is a solution to the Hamilton-Jacobi equation. It is argued that this wave function 
corresponds to an ensemble of classical trajectories satisfying the first integral p = V-S, 
with |Cp giving a measure on the ensemble. Although probably correct it is somewhat 
heuristic and can only be used in the oscillatory regime. (See Ref.[7], for example, for a 
discussion of these issues). 

Recent more successful work with the Wheeler-DeWitt or similar equations involves the 
induced inner product (also known as Rieffel induction or refined algebraic quantization) 
[8,9]. Here, one considers eigenvalues of the Wheeler-DeWitt operator, 

H\^Ek)=mE'k') (1-5) 

where k is the degeneracy label for each E. The spectrum is typically continuous in E in 
which case the states are normalized via Eq.(1.4) according to 

{<i^Ek\'^E'k') = S{E-E')6kk' (1.6) 
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and one can now see why (^'l^') diverges. The induced inner product is then, loosely- 
speaking, to drop the factor 5{E — E') as E and E' are set to zero. This procedure can 
be defined rigorously and induces an inner product on the zero energy eigenstates. (This 
procedure is not of course necessary when the spectrum is discrete). 

Related to these issues is the prevalent idea that any operations performed on the wave 
function in the computation of physically interesting probabilities should commute with H 
[10,3,11]. Mathematically, this is to respect the symmetry of the theory, reparametrization 
invariance, expressed by the constraint equation (1.1). Physically, it is connected with 
the fact that the universe is a genuinely closed system, and all realistic measurements are 
carried out from the inside, so cannot displace the system from a zero energy eigenstate 
of H. 

Given these preliminaries, turn now to the questions one would like to ask of the wave 
function of the system in order to extract useful cosmo logical predictions from it. We 
are interested in the notion that the wave function corresponds in some way to a set of 
trajectories. Let us therefore ask the question. 

What is the probability that the system is found in a series of regions in configuration 
space, Ai, A2,---An? 

Note that the question is stated in such a way that does not involve time. There is no 
requirement that the system enters one of the regions at a particular "time", or that 
the regions are entered in a particular order. We cannot ask this because the Wheeler- 
DeWitt equation does not know about such an ordering parameter. In the classical case 
the corresponding situation consists of a statistical ensemble of classical trajectories with 
the same fixed energy. The trajectories are simply curves in configuration space and it 
is straightforward to determine the probability that a given curve passes through a given 
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region at any stage in its entire history. The question is more involved in quantum theory, 
since quantum theory is somewhat resistant to the notion of a trajectory (in the non- 
relativistic case, it involves specifying positions at different time, which do no commute). It 
is, nevertheless, important to develop this notion, since the timeless nature of the Wheeler- 
DeWitt equation cries out for an interpretation in terms of entire histories of the universe. 
The aim of this paper is to offer one possible way of giving meaning to the above question 
in the quantum theory of simple cosmological models. 

Intuitively, one would expect that the question can be formulated and answered using a 
simple toolbox of parts: the quantum state satisfying the constraint, projection opera- 
tors onto the regions A^^, or maybe projections onto some class of operators which commute 
with the constraint H. One might also expect to find the Green function associated with 
the Wheeler-DeWitt equation, which has the form, 

roo 1 

G(x,y) = i dr (x|e-^(^-^^)ny) ^ (x|— — -|y) (1.7) 
Jo {H - ^e) 

(this is the analogue of the Feynman propagator) . It might possibly also involve one of the 
other types of propagators obtained by integrating r over an infinite range in this expression 
[12,13,14]. The question is then exactly how one stitches all these components together to 
make a plausible probability distribution describing trajectories passing through a series 
of regions. 

The decoherent histories approach offers an approach to answering this question and it 
does indeed use some of the above elements [15]. It is particularly adapted to this sort of 
situation since it directly addresses the issue of defining a quantum notion of "trajectory" 
and this approach is currently being investigated in this context [16] (see also Ref.[l]). 
Other approaches involving observables - operators commuting with the constraint - have 
also been considered [17,18,9,19]. Most importantly, Kiefer and Zeh [20] and Barbour 
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[21,22,4] have devoted much effort to elucidating the emergence of trajectories and of time 
from the timeless Wheeler-DeWitt equation. 

The approach we adopt here stems from Barbour's observation [22] that a substantial 
insight into the Wheeler-DeWitt equation may be found in Mott's 1929 analysis of alpha- 
particle tracks in a Wilson cloud chamber [23]. Mott's paper concerned the question of 
how the alpha-particle's outgoing spherical wave state, e^^^/R, could lead to straight line 
tracks in a cloud chamber. His explanation was to model the cloud chamber as a collection 
of atoms that may be ionized by the passage of the alpha-particle. They therefore act 
as measuring devices that measure the alpha-particle's trajectory. The probability that 
certain atoms are ionized is indeed found to be strongly peaked when the atoms lie along 
a straight line through the point of origin of the alpha-particle. 

Although Mott seems to have had in mind a time-evolving process, he actually solved 
the time-independent equation 

{HQ + Hd + \Hint)m = E\^i>) (1.8) 

Here is the alpha-particle Hamiltonian, Hfi is the Hamiltonian for the ionizing atoms, 
and Hi^i describes the Coulomb interaction between the alpha-particle and the ionizing 
atoms (where A is a small coupling constant). Now the interesting point, as Barbour 
notes, is that Mott derived all the physics from this equation with little reference to time. 
Mott's calculation is therefore an excellent model for many aspects of the Wheeler-DeWitt 
equation. Barbour has elucidated this very eloquently, showing how it sheds light on a 
number of different aspects [22,21,4,24]. 

Barbour's discussion is largely qualitative. The point of the present paper, by contrast, 
is to extract quantitative information from the comparison between the Mott calculation 
and the Wheeler-DeWitt equation. Mott derived the straight line tracks by looking at the 
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wave function associated with two atoms being in the ionized state, for the special case of an 
outgoing wave initial state. But since this is elementary quantum mechanics, it is a simple 
matter to generalize it to arbitrary initial states and other types of detectors models, and to 
derive a detailed expression for the probability distribution. Therefore, Mott's calculation 
points the way towards a general expression for the probability distribution for the system 
passing through a series of regions in configuration space without reference to time. This 
is what we will work out in detail in this paper. 

We consider a system in an N dimensional configuration space with coordinates x 
described by a Hamiltonan Hq, which may be of the form (1.2), but the simplest case we 
consider is a free particle. It is coupled to a set of detectors via an interaction H^nt and the 
state of the whole system is given by the solution to (1.8). Mott used the electronic degrees 
of freedom of atoms as detectors. However, the essence of the calculation is maintained 
with a much simpler detector model. The detector we use consists of a two state system 
with H(i = and detector states |0) and |1), where, 

a|0) = 0, a|l) = |0), at|0) = |l), at|l) = (1.9) 

We take 

Hint^Yl //tW(«fc + 4) (1-10) 
k 

Here, /fc(x) is spatially localized in the region A^^. One could, for example, take fj^ 
to be a window function which is 1 in Aj^ and outside it, but we will not restrict 
to this choice. If the detector in A^^ is "initially" in the ground state |0), it will be 
displaced into the excited state |1) if the particle's trajectory x(t) enters A^^ and stays 
in the ground state otherwise. Of course, in the timeless context of the Wheeler-DeWitt 
equation, "initially" has no meaning. Instead, following Mott, the appropriate condition to 
impose is that the detector is in the ground state in the absence of coupling to the system. 
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This detector is far from realistic, not least of all because it can return to its ground state 
if the particle spends too much time in the detector region. We will discuss its problems 
and possible improvements. We note however, that similarly simple detector models have 
been profitably used elsewhere, e.^., in the Coleman- Hepp model [25]. (See also Ref.[26]). 

It is perhaps worth noting that the question considered here bears a close resemblance 
to the arrival time and tunneling time questions in non-relativistic quantum mechanics 
[27,28]. There also, time enters in a non-trivial way, and equivalent classical approaches 
to the problem are inequivalent at the quantum level. A variety of approaches have been 
brought to bear on these problems, including, as here, explicit detector models. 

In Section 2, we solve the system Eq.(1.8) using the simple two state detector model. 

Using the results, we then ask, in Section 3, some simple questions of the detected 
wave function. Does the detected wave function still obey the Wheeler-DeWitt equation? 
We find that it does exactly outside the detector region, and that it does approximately 
(for small A) inside the detector region. We compute the probabilities for detection and 
see, as Mott essentially saw, that they are strongly peaked when the detectors lie along 
a classical trajectory. We also observe that the resulting amplitude for detection bears a 
very close resemblance to a formula written down by Hartle [29] (without using an explicit 
detector model) and we discuss the connections with his result. We also discuss various 
other aspects of the solution in relation to detection and the timelessness of the solutions. 

The results of Section 3 indicate that the wave function of the system may, in some 
sense be regarded as a superposition of states each of which is peaked about an entire 
classical history. To demonstrate this explicitly, we specialize, in Section 4, to the case of a 
collection of harmonic oscillators and introduce a new type of coherent state, the "timeless 
coherent states" . These states are eigenstates of the Hamiltonian, and are therefore time- 
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independent, but are peaked about classical trajectories. Any eigenstate may be expanded 
in terms of theses states, and we show that a series of detections along a classical path 
essentially projects the state down onto a timeless coherent state. 

Since the detector is so simple, its dynamics may be solved exactly and this is carried 
out in Section 5. This calculation confirms that the detector model is only physical realistic 
in the perturbative regime (when the particle spends only a short time in each detector 
region). Although the solution is exact and leaves the boundary conditions general, it turns 
out that it is not very useful for the Mott solution, since the boundary conditions lead to a 
rather inelegant solution. On the other hand, it is by no means clear that one is required 
to take the Mott boundary conditions for the analagous situation in quantum cosmology, 
and given the freedom to choose different conditions, an elegant alternative solution for 
the detector amplitude is obtained. It is in fact almost the same as the amplitude Hartle 
proposed [29]. 

In Section 6, we briefly describe a more elaborate detector model, in which the detector 
is a simple harmonic oscillator coupled to the particle with the same coupling (1.10). The 
solution has a nice path integral representation and suffers fewer shortcomings that the 
two-state detector. It also clearly illustrate the peak about classical paths. 

We summarize and conclude in Section 7. 

2. DETECTION AMPLITUDE FROM 
THE TWO-STATE DETECTOR 

We now solve the eigenvalue equation 
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(2.1) 



with H given by Eq.(1.8), and the detector is the simple two state detector described in 
Section 1, with = 0. We will for convenience refer to this equation as the Wheeler- 
DeWitt equation (and it is convenient to retain a non-zero value of E). It now reads 

{H - E)m = {Ho + XH.^t - Em = (2.2) 

We initially use only two detectors, so 

2 

Hint=J2 /fcW(«fc + 4) (2-3) 
k=l 

but the generalization to a arbitrary number of detectors is straightforward. We solve 
perturbatively by writing 

l^r) = l^r(O)) + A|*W) + X^l^C^)) + ... (2.4) 

We require that in the absence of coupling to the detectors, the detectors are in the state 
of no detection, |0). This means that 

|*(0)) = 1^^)10)10) (2.5) 

and I';/') is the state we are trying to measure. Inserting in (2.2) and equating powers of A, 
we get, 

{Ho - E)\^^^^) = (2.6) 
{Ho-E)\^(^^)^-Hm\^(^^) (2.7) 
{Ho-E)\^('')) = -Hi^t\^(')) (2.8) 

and similarly to higher orders. The first relation says that lip) must obey the unperturbed 
eigenvalue equation, as expected. Inserting (2.5) into Eq.(2.7), we get 

{Ho - i?)|^'W) = -M^mim - /2(x)|V')|0)|l) (2.9) 
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This is readily solved by writing 

1*^'^) = l*S)|0)|0) + |*5i^)|l)|0) + 1*^1^)10)11) + |*5ft|l)|l) (2.10) 



and we discover that 



{Ho 


-E)|*, 


{Ho 




{Ho 


-E)|*, 


{Ho 


-E)\^ 



00 / 

lil / - 

(IX _ 

01 / - 



(2.11) 
(2.12) 
(2.13) 
(2.14) 



Eqs.(2.12), (2.13) may be solved with the assistance of the Green function G, defined by 
Eq.(1.7) (with H replaced by Hq — E). It obeys the equation 

{Hq - E)G = 1 (2.15) 

(For convenience we use an operator notation in which G is the operator with coordinate 
representation G(x, y) = (x|G|y) and the right-hand side of (2.15) would be the delta- 
function (5(x, y) in the coordinate representation). We then obtain 

l*S) = -G'/i(x)|V^) + |0l) (2.16) 
|*£^) = -G/2(x)|V^) + |02) (2.17) 

where |0i,2) ^ire solutions to the homogeneous equation 

(i/o-^)l01,2) =0 (2.18) 

To fix some of these solutions more precisely, we need to appeal to boundary conditions. 
This is a subtle issue and depends very much on the precise context. Mott was concerned 
with the particular case of an outgoing spherical wave and imposed boundary conditions 
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appropriate to this case. This led him to set |</>i,2) and l^n ) to zero (since otherwise it 
represents a stream of incoming particles fired at an already excited detector) [23]. We 
are not obviously compelled to make the same choice of boundary conditions in quantum 
cosmology, and we will return to a discussion of this important issue in Section 3(C). But 
for the moment, we work with the Mott solution. 

At first order only one detector is stimulated into the excited state. The system wave 
function correlated with detector state |1) now is 

|V'l) = -AG/i|V) (2.19) 
and the probability that the detector registers is therefore 

PI = (V'llV'l) = X^miG^GfiliP) (2.20) 

(When the spectrum of Hq is continuous, expressions of this type need to be regularized 
along the lines of Eq.(1.6), but we will carry out this explicitly only when we need to 
calculate it in more detail). 

To get two detectors to register, we need to go to second order. We now have, from 
(2.8), and the solution (2.10), 

{Hq - E)|*(2)) = _ (^f^^ai + a\) + /2(a2 + 4)) 

= /2G/1 1^^)11)11) + hGf2\mm + ■■■ (2.21) 

where the omitted terms on the right-hand side are proportional to |0)|0), |0)|1) and |1)|0), 
and will not be needed. Again we may solve by expanding as in (2.10). Here we write 
down only the term required, which describes two detectors being excited 

|^r(2)) ^ + (2.22) 
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and it is readily seen that the solution is 

1*1?) = {Gf2Gfi + G/iG/a) IV) (2.23) 

Again, following the spirit of the Mott solution, possible homogeneous solutions are set to 
zero. We now have the system wave function correlated with two detectors registering: it 
is 

M = a2 {Gf2Gfi + GfiGf2) IV) (2.24) 
and the probability is (V2|'02)- 

The analysis is readily extended to an arbitrary number of detectors, but it is easy to 
anticipate the result from (2.24): for n detectors, the amplitude is 

IVn) = A" {GfnGfn-i ■ ■ ■ G/2G/1) |V) + symmetrizations (2.25) 

where "symmetrizations" means add all possible permutations of 1, 2,3, ■■ - n. It is clear 
that these terms are there to ensure that there is no preference in the order in which 
each of the detectors registers, reflecting the genuinely timeless nature of the underlying 
dynamics. Eq.(2.25) is the main result of this Section. 



3. PROPERTIES OF THE SOLUTION 

We have shown that the wave function for the whole system when there are, for exam- 
ple, two detectors, takes the form 

I*) = |Vo)|0)|0) + iVi) (|1)|0) + |0)|1)) + |V2)|1)|1) (3.1) 

We can now ask various questions of the detected wave function IV2)) or more generally, 
Eq.(2.25). 
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3(A). Does the Detected Wave Function 
obey the Wheeler-DeWitt Equation? 



As stated in the Introduction, a prevalent idea in quantum gravity is that all observables 
should commute with the total Hamiltonian [10,3,11]. Related to this is the notion that 
"measurements" of the wave function (whatever this may mean in general) should not 
displace the system from its eigenstate of the Hamiltonian. Given that we have presented 
here an explicit model of detection, it is perhaps of interest to ask to what extent this idea 
holds up. 

We have, using (2.15), and taking the simple case n — 2, 

{Ho - E)\^^) = a2 (/2G/1 + /1G/2) IV') (3.2) 

In configuration space, the right-hand side zero, except in the detector regions Ai, A2, 
because the functions fi and /2 are localized there. In these regions, it is of order A^, which 
we regard as small in comparison to the terms on the left. Hence the measured wave func- 
tion approximately obeys the Wheeler-DeWitt equation. This is no surprise. In standard 
quantum mechanics, a physically measured system does not obey the Schrodinger equa- 
tion but has corrections due to the measuring device. That it obeys the Wheeler-DeWitt 
equation only approximately is not in conflict with exact reparametrization invariance, 
since the wave function for the entire system always obeys the Wheeler-DeWitt equation 
exactly. On the other hand, one wonders whether it might not be possible to devise a 
detection scheme in which the detected amplitude obeys the Wheeler-DeWitt equation 
exactly. For example, given the simple toolbox of parts outlined in the Introduction (such 
as the Green function (1.7)), it would not be unreasonable to guess that the detection 
amplitude might be of the form (2.25) but with the G given by Eq.(1.7) replaced by the 
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one obtained by integrating r over an infinite range. This gives a solution to the Wheeler- 
DeWitt equation, and as a consequence, the modified detection amphtude would also obey 
the Wheeler-DeWitt equation exactly. 

3(B). Amplitudes for Classical Paths 

Of greater interest is the question of the configurations about which the amplitude 
(2.25) (or the associated probability) is peaked. The amplitude may be written in an 
explicit coordinate representation as 

(x^lV^^) = A^y" d^x„---d^X2d^xi G(xy,x„)/„(x„) 

X G(xn,x„_i)/„_i(x„_l) ••• G(x2,xi)/i(xi) V(xi) (3.3) 

(plus symmetrizations) . It has the form of approximate projections onto the regions Aj^ 
(exact projections if the ff^ are window functions) with evolution between regions described 
by the fixed energy propagator G(xy^_|_x, x^). It is analogous to the amplitude for a history 
of positions at a sequence of times in non-relativistic quantum mechanics, which is known 
to be peaked about classical trajectories [30] . But note that here the evolution is at fixed 
values of the energy, and there is no reference to time. 

We can estimate the form of (3.3) using a WKB approximation for the fixed energy 
propagator [31]. It is given by an expression of the form, 

G'(x",x') = C(x",x')e^^('^"''^') (3.4) 

where C is a slowly varying prefactor and <S'(x",x') is the fixed energy Hamilton-Jacobi 
function, i.e., the action of the classical solution from x' to x" with fixed energy E [32]. 
The initial and final momenta of this classical solution are 

p" = Vx"-S(x", x'), p' = -Vx'^(x", xO (3.5) 
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In terms of (3.4), the amphtude (3.3) may be written, 

/n 
d^Xn • • -d^xsd^xi Yl C(xfc+i,Xfc)/fc(xfc) 
k=l 

X exp ^^(xyt+i,Xfc)j '(/'(xi) (3.6) 

where x„_|_i = xj. 

Consider the integrals over X2 • • • x,^ with xi (and x„+i) fixed. Suppose for the moment 
that the functions fj^ are absent so the integrals are unrestricted. By the stationary phase 
approximation, the dominant contribution to the integral comes from the values of X2 • • • x^ 
for which the phase is stationary, i.e., for which 

n 

Vx,$^5(x^+l,xfc) = (3.7) 
k=l 

ior j — 2, ■ ■ ■ n. This equation means that 

Vx,-S(xj+i, x_,-) + Vx,5(xj, xj_i) = (3.8) 

Using Eq.(3.5), this implies that the final momentum of the classical path from ^j-i to Xj 
is equal to the initial momentum of the classical path from Xj to ^j+i- It is not difficult 
to see that this in turn implies that the point Xj must lie on the classical path from y^j-i 
to Xj+i- Hence, the stationary phase points of the whole integral (3.6) lie on the classical 
path from xi to x^. The approximate value of the integral is the integrand of (3.6) with 
the stationary phase point values inserted. 

Now consider what happens when the restricting functions are present. If the regions 
(about which the f]^ are concentrated) include the stationary phase points (and if the 
regions are larger than the fluctuations about these points), then, since the integral takes 
it dominant contribution from these points, the presence of the /j^ makes little difference 
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and the integral is given once again by its stationary phase value. On the other hand, if 
one or more of the fj. lie far away from the stationary phase points, then, since the integral 
is prevented from taking a contribution from these points, its value must be much smaller 
than the stationary phase value. We thus see that the amplitude (3.3) will be largest when 
the regions Aj^ are chosen to include the stationary phase points of the integral. As we 
have shown, the stationary phase points lie along the classical path from xi to xy. It 
follows that the amplitude (3.2) will be largest when the regions Af^ are chosen to lie along 
a classical path. 

Mott's argument for straightline paths in an expression analogous to (2.25) relied on 
the explicit from of the Green function in the free particle case, and on a special initial state 
[23] . Here we see that the peaking about classical paths can be seen, at least heuristically, 
from elementary properties of path integrals, for a broad class of Hamiltonians and initial 
states. Bell has also discussed the Mott calculation at some length [33]. He notes that 
the first projector fi spatially localizes Mott's initial wave function, but in a realistic 
ionizing event, the resultant uncertainty in momentum can still be extremely small. As 
a consequence the angular spread of the wave packet in its subsequent evolution can be 
extremely small, hence the appearance of a straight line track. 

3(C). Comparison with a Result of Hartle 

The result (2.24) is very closely related to a result of Hartle [29]. He considered a 
simple model quantum cosmology with a Hamiltonian quadratic in the momenta, as here, 
and asked for the amplitude that the system passes through two regions of configuration 
space, Ai, A2. Using some simple arguments about propagators and elementary principles 
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of quantum theory, he showed that the amphtude is (in the language of the present paper) 

(G/2G/1 + G/iG/s) IV;) - (Gt/2Gt/i + Gt/iGt/2) 1^) (3.8) 

where ip is a solution to the Wheeler-DeWitt equation and fi and /2 are taken to be exact 
projectors onto the regions Ai, A2. Other than the factor of in (2.24) which is not 
important, Hartle's result differs from (2.24) by the subtraction of an identical term but 
with G replaced by G^ (which is generally not the same as G). Hartle argues that this 
should be there on the grounds that, in an expression like (2.25) with G represented by 
(1.7), the time parametrization should not have a preferred direction with y the "initial" 
point of the parametrization and x the "final" point, hence we should sum the amplitude 
over both possible parametrization directions. 

To understand why this term can be there in the present calculation, let us first say 
a little more about the Green function. G, as defined by Eq.(1.7), may be written more 
explicitly as 

G(x,y,£;) = y (3.9) 

^ ' ^ E - En + ie ^ ^ 

n 

where ttn(x) are eigenfunctions of Hq with eigenvalue En- G{x,y,E) is real when E lies 
in the discrete part of the spectrum and complex when E lies in the continuous part (see, 
for example, Ref.[34]). Hence in the free particle case considered by Mott, G ^ g\ , but for 
the harmonic oscillator (considered in Section 4) G = G^ . Quantum cosmological models 
usually have a spectrum which is at least in part continuous, so we expect G ^ G^ 'm. 
general. 

Now recall the comments after Eq.(2.18), where it was noted that we are by no means 
obliged in quantum cosmology to take the same boundary conditions as Mott. Since 
also satisfies Eq.(2.15), we may use it in place of G to generate solutions to the detector 
amplitude. Because G — G^ = i5{HQ — -E"), it is easily seen that the difference between 
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using G and is a homogeneous solution in (2.16), (2.17). It is also true that, in higher 
order perturbations, we may use G or or some combination. Therefore, there is a more 
general class of solutions for the detection amplitude which are sums of terms of the form 
(2.25) with some of the G's replaced by G^^s (and with a suitable overall normalization). 
In particular, Hartle's amplitude falls into this enlarged class of solutions, so there is no 
conflict with his result. 

Mott made a particular choice of solution appropriate to the physical situation he 
as investigating. In the case of relativistic field theory in Minkowski space, one would 
normally impose some sort of causality requirement to fix the solution more precisely, and 
so to choose between G and (since G is essentially the Feynman propagator). One 
could require, for example, that the wave function for the whole system is affected by 
the detector only in the future light cone of the detector region. In quantum cosmology, 
however, although the metric has hyperbolic signature like Minkowki space, it is by no 
means clear that one is obliged to impose an analogous requirement, and in fact it is 
difficult to see exactly how to do this in general since the configuration space is usually 
not globally hyperbolic. (See, however, Ref.[12]). Instead, one might expect to fix the 
choice of G or G^ by appealing to cosmological boundary conditions. The no-boundary 
proposal of Hartle and Hawking, for example, picks out a wave function that is real [35] . 
It is sometimes claimed that this corresponds to a "time-symmetric" wave function [36]. 
It therefore stays most closely to the timeless nature of the Wheeler-DeWitt equation and 
in some sense represents a complete abandonment of any fundamental notion of causality. 
It is now interesting to note that Hartle's detection wave function (3.8) is in fact real if 
i/j is real. In many ways this therefore seems like the most natural solution to take in the 
case of the Wheeler-DeWitt equation. We do not, however, in this paper commit to any 
particular choice of boundary condition. 
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In summary, therefore, a more general solution to the detector model is a sum of terms 
of the form (2.25) involving both G and , and this more general solution includes Hartle's 
result (3.14). Note also that the replacement of G by does not affect the discussion of 
the peaking of the amplitude about classical trajectories. 

3(D). On Timelessness and Detection 

It is perhaps worth elaborating on a feature of the detector model which appears at 
first sight to be incompatible with the timelessness of the Wheeler-DeWitt equation. We 
have coupled the system to a series of detectors via the interaction Hamiltonian (1.8). 
This Hamiltonian describes a situation in which, along a trajectory x(t), the detector is in 
the ground state "before" the trajectory enters the detector region and in an excited state 
"after" it has passed through the region. Along a classical trajectory x(t) in which there 
is a notion of time, and of before and after, this is undeniably correct. (The parameter t 
simply labels the points along the curve - it may be take to be, for example, the distance 
along the curve from some reference point.) But how are we to understand how the detector 
works in the genuinely timeless world described by the Wheeler-DeWitt equation? There 
is no before and after and there is no preferred direction of time. 

The above results effectively imply that each solution to the Wheeler-DeWitt equation 
may be regarded as a superposition of states each of which is concentrated along an entire 
classical trajectory in configuration space (and we will see this in more detail in Section 
4). What then seems to be happening in the present model including the detector, is that 
each trajectory carries a label indicating whether or not it passes through the detector 
region at any stage along its entire length. In the (perhaps restrictive) language of time, 
at each point along the trajectory, the label allows the trajectory to "know" whether it 
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passed through the detector region in the past, or will pass through it in the future. But 
it is only on adopting this temporal language that the situation seems paradoxical. The 
paradox vanishes when ones speaks the vocabulary of entire trajectories in configuration 
space, and one can see this in the solution (3.1). The wave function of the entire system 
is written as a correlated state in which the state correlated with the detectors in the |1) 
state is the state (2.25): the detected state is indeed concentrated on trajectories that pass 
through Ai, • • • 

The issues discussed in this Section may be of relevance to the perennial debate on the 
question of time asymmetry in quantum cosmology [37,20]. 



Since the peaking about classical paths is the most important property of the amplitude 
(2.25), it is worth exploring it in more detail for the special case of a collection of harmonic 
oscillators, where it is possible to show very clearly how the solutions to the Wheeler- 
DeWitt equation correspond to superpositions of states peaked about classical paths. We 
will introduce a class of coherent states appropriate to the timeless theories considered here 
and which are natural analogues of the standard coherent states of the harmonic oscillator. 

The Hamiltonian for a set of N identical harmonic oscillators is 



4. COHERENT STATES FOR TIMELESS DYNAMICS 




(4.1) 



In this case the spectrum of Hq is discrete and we have 




(4.2) 
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Since S{Hq — E) is now in fact a true projection operator we may write 5{Hq — Ey = 
S{Hq — E) without having to worry about regularization through the induced inner product, 
as in the continuous case. (In this expression E is allowed to take only the discrete values 
corresponding to the spectrum of Hq). The Green function G is given as before by (1.7) 
with H replaced by Hq — E. For a one- dimensional oscillator, Eq.(4.2) is equivalent to 



where l^') is the energy eigenstate. In more than one dimension the energy eigenstates are 
degenerate, so Eq.(4.2) has the form 



d 

where \E, d) are the energy eigenstates with degeneracy label d. 

The standard coherent states (see Ref.[38], for example) are denoted |p,x) and they 
have the important property that they are preserved in form under unitary evolution. 



where pt,Xi are the classical solutions matching p, x at t = 0, hence they are strongly 
peaked about the classical path. We are interested in finding a set of states which are 
analogues of these states for the timeless case. That is, they should be eigenstates of Hq^ 
and should be peaked about the classical paths of given fixed energy in phase space. It is 
not difficult to see that a set of states doing the job are. 



5{Hq-E)^\E){E\ 



(4.3) 




(4.4) 




(4.5) 



|(/>px)=5(ifo-^)|p,x) 




(4.6) 
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These states are not in fact normalized to unity but we shall see that it is useful to work 
with them as they are. Since the states |pi,Xi) are concentrated at a phase space point 
for each t, clearly integrating t over a whole period produces a state which is concentrated 
along the entire classical trajectory. Each state is labeled by a fiducial phase space point 
p, X which determines the classical trajectory the state is peaked about. Under evolution 
of the fiducial point p, x to another point, Ps, Xg, say, along the same classical trajectory, 
the state changes by a phase, 

|0px) - = e^^^|0px) (4.7) 

as may be seen from (4.5), (4.6). We will refer to these states as timeless coherent states. 
Their properties are in fact very similar to the usual coherent states. 

Note that the coherent states |p, x) are in fact already approximate eigenstates of Hq, 
with eigenvalue ^(P^ + x^); ^ long as |p|, |x| are much larger than the coherent state 
quantum fiuctuations. Given E, it therefore seems reasonable to choose the values p, x in 
the fiducial coherent state so that E = ^(P^ + x^)) when constructing the timeless states 
(4.6). 

Two timeless coherent states of different energy are exactly orthogonal. The more in- 
teresting case is that in which they have the same energy, and then they are approximately 
orthogonal if they correspond to sufficiently distinct classical solutions. This is because we 
have 

(0p'x'|0px) = {p',^'\S{Ho-E)\p,^) 

= I £e^^^(p',x'|pi,x,) (4.8) 

From the properties of the standard coherent states we know that 

|(p',x'|pt,xt)| <1 (4.9) 
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with equality if and only if p' = pt and x' = x^. Moreover, the overlap of two coherent 
states is exponentially small if they are centered around phase space points that are suffi- 
ciently far apart. It follows that if p',x' does not lie on, or close to, the trajectory pt,^ti 
the overlap (p', x'jp^, x^) will always exponentially small for all t. The integral over t in 
(4.8) will then give a result that is much smaller than the case in which p', x' does lie 
on, or close to, the trajectory pt,:x-t (because in the latter case the overlap (p', x'|pf, x^) 
becomes close to unity for some value of t). The timeless coherent states are therefore 
approximately orthogonal for sufficiently distinct classical trajectories. 

The standard completeness relation for the coherent states is 



Multiplying both sides by d{HQ — E) from the left and right, and using (4.6), we get 



Since 5{Hq—E)\iI)) = on any solution to the eigenvalue equation [Hq—E)\iI)) = 0, this is 
as good as a completeness relation on the set of solutions to the eigenvalue equation (which 
is all we are interested in). We may therefore write any solution |?/^) as a superposition of 
timeless coherent states. 



It is then tempting to interpret (^pxlV') ^ ^be amplitude that a system in the state 
will be found on the classical trajectory corresponding to the timeless coherent state |(/>px)- 
Again using the fact that 5{Hq — E)\'^) = lip) it is easy to see that this is in fact the same 
as (p,x|7/;), which is the amplitude for finding the system at the phase space point p,x 
labeling the trajectory. The probability is then simply |(p,x|'0)p. This is a simple and 




(4.10) 




(4.11) 




(4.12) 
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intuitively appealing result: the classical trajectory is completely fixed by its initial values 
X, p, hence we expect that the probability for being found on a certain classical trajectory 
is the same as the probability for being found at the initial phase space point that labels 
it. 

This interpretation is put forward with a small note of caution, however, since the 
sum over p, x in Eq.(4.12) is not only over states which are only approximately orthogonal 
(as with the usual coherent states), but, because of the property (4.7), includes some 
redundancy in the summation. In particular, \(f)p-x) {(l>px.\ is invariant along the classical 
phase space trajectory of its fiducial point. Since this only produces some sort of constant 
factor, it may not make any difference, and indeed, the above interpretation appears to 
produce sensible results. Still, it would be desirable to include, if possible, some sort of 
"gauge fixing" which factors out this redundancy. This will be explored elsewhere. 

Given all of this background, we may now reconsider the detector amplitude (2.25). 
We expand the initial state IV') in the timeless coherent states, as in Eq.(4.12). When fi 
operates on the state, it has the effect of restricting the sum over timeless coherent states 
to only those that pass through the region Ai, but otherwise leaves the timeless coherent 
states preserved in form (as long as the region Ai is larger than the width of the coherent 
states). Next, to consider the operation of G, we write it as 



where H = Hq—E. When the timeless coherent states are operated on with the exponential 
factor in (4.13) it has essentially no effect, since H\(f)px) = 0. (The time integration leads to 
a constant factor, which may, however, need regulating along the lines of (1.6)). Operating 
with /2 further restricts the sum to those states that pass through both Ai and A2, and 
then subsequent operation G again does essentially nothing. And likewise to the end of the 




(4.13) 
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chain. We therefore find that the detector amphtude (2.25) consists of a superposition of 
only those timeless coherent states that pass through all the detector regions Ai, A2 ■ ■ ■ A„. 
This could be quite a large sum of states if the regions Aj^ are large. However, it will consist 
of essentially just one timeless coherent state if the detector regions lie along a classical 
trajectory and if their size is just bigger than the spatial width of the wave packet. The 
detection amplitude along this classical trajectory is then equal to (^pxlV')) agreement 
with the analysis based on Eq.(4.12). 

Rovelli has written down a coherent state of the type considered here [17] in the 
context of a very similar model, although its properties were not explored and exploited 
as they are here. Klauder, in his approach to the quantization of constrained systems 
using coherent states, considered the projection of the standard coherent states onto the 
contraint subspace, hence in essence wrote down states of the form (4.6) [39]. He did 
not, however, consider their use as an interpretational tool. Wave packet solutions to 
the Wheeler-DeWitt equation which aproximately track the classical trajectories for more 
interesting cosmological models have been considered in Refs. [40,41]. 



5. EXACT SOLUTION TO THE DETECTOR DYNAMICS 



In this Section we discuss an alternative method of solution of the eigenvalue equation 
(2.2). As we shall see, it does not in fact give a very elegant representation of the Mott 
solution, which is why it was not used above. If, however, one is permitted to use different 
boundary conditions, as may be reasonable in quantum cosmology, then it provides an 
alternative possible solution to the detector dynamics. 

We consider first the case of a single detector. The key observation is that solution to 
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the eigenvalue equation (2.2) may be generated using either of the expressions, 



I*) = 5{H-E)\(j)) 
~ 2^ 



^-iT{Ho+\Hint-E)\ 



(5.1) 

(5.2) 



oo 



Here |0) is an arbitrary fiducial state in the joint system-detector Hilbert space. It is 
ambiguous up to the addition of a term of the form [H — E)\(j)'). For the moment we keep 
10) general and take 

!</>) = |X0)|0) + |X1)|1) (5.3) 
To evaluate (5.1) or (5.2) we introduce the eigenstates of a + at, which are 

1 



V2 



(|0)±|1)) 



(5.4) 



The expression (5.1) is then readily evaluated with the result 



l*> = lV'nrf)|0) + |V'd)|l> 



(5.5) 



where. 



\^nd) = \ m + A/i) + S{h - A/i)) Ixo) + I m + A/i) - 6{h - A/i)) Ixi) (5.6) 
M = I iS{h + A/i) - 6{h - A/i)) Ixo) + I {S{h + A/i) + 6{h - A/i)) |xi) (5.7) 

(recalling that h = H — E). These expressions are perhaps more easily appreciated using 
a path integral respresentation: 

JtE 



(5.8) 



^„rf(x;) = / dr e*^^ / V^{t) exp (^5[x(t)]) 



— oo 



X 



cos^A^ (^^/l(x(0)^ Xo(xo) +«sin ^A^ (^^/l(x(t))^ Xl(xo) 



/oo r 
dre^^^ / Px(t)exp(z5[x(t)]) 
-00 J 



X 



isin (xj\th{^{t))^ xo(xo) +COS (^X dth{yi{t)) ) xi(xo) 
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(5.9) 



Eqs.(5.8), (5.9) represent the exact solution to Eq.(2.2) with one detector in place. 

Turn now to the question of the fiducial state. The condition (2.5), suggests that we 
should take 



where \'^) = S{Hq — E)\xo), and hence that |xi) = 0. In fact, in the induced inner product 
scheme, we may take |xo) = I'*/') 7 since effectively 6{Hq - Ef = 5{Hq - E). This therefore 
yields a path integral expression for the amplitude for detection and has the property that 
the factoring condition (2.5) is satisfied when A = 0. So let us first explore the properties 
of (5.8), (5.9) with xi = 0. 

The nature of the sums over paths x(t) in Eq(5.8) and (5.9) is governed by the quantity 



appearing in the sine and cosine factors. With fi normalized to be dimensionless, this 
quantity has the dimensions of time, and is essentially the time spent by the path x(t) in 
the region Ai around the detector. (This is not of course a physically measurable time. 
The paths x(t) in the path integral have a well-defined notion of time, but after summing 
over the total time duration of each path r the final result is time-independent). Now, in 
the amplitude for no detection, (5.8), the factor cos{Xt^) is 1 for = 0, and decreasing for 
increasing from zero. In the path integral, it therefore has the effect of suppressing paths 
that pass close to the detector, and favours paths that stay away from it. Similarly, in the 
detection amplitude, (5.9), the factor sin(AT(^) is zero for paths that spend no time near the 
detector, and non-zero for paths that enter the detector region. It therefore enhances the 
amplitude for paths entering the detector region. The sine and cosine factors therefore, in a 
very crude way, enforce restrictions on the paths corresponding to entering or not entering 
the detector region, as one would expect. We can also see, however, that these factors only 



= lxo)|o) 



(5.10) 




(5.11) 
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do their job well if Ar^^ is somewhat smaller than 1, indicating that the detector model is 
only physically sensible in the perturbative regime about A = 0, as expected. 

It is now important to check the agreement between the exact result above and the 
perturbative result of Section 2. It is not difficult to show that the exact solution (5.9) 
with xi = Q does not in fact agree with the perturbative solution (2.19). It differs by the 
presence of homogeneous solutions in the small A limit of (5.9). We will not go into details 
but it may be shown using (5.7) and the identity, 

S{h + A/) = 6{h) + iX (cfG - Gt/ct^ + 0(a2) (5.12) 

(which is proved using the exponential representation (5.2)). Recall that in the pertur- 
bative solution the homogeneous solutions were removed at each order in perturbation 
theory essentially by inspection. Since no corresponding condition has been imposed here, 
it is not surprising that these spurious solutions crop up. We can, however, get agreement 
if we choose |xi) = ■iA/G^|xo) in Eq.(5.9) (again proved using (5.12), although there is 
no obvious independent reason for making this choice. Furthermore, the presence of the 
Green function in the fiducial state rather destroys the elegance of the path integral 
representation compared to the case xi = 0- 

On the other hand, although the Mott solution is not readily obtained, the removal 
of the homogeneous solutions in the perturbative solution is a subtle matter of boundary 
conditions, above and beyond the basic factoring condition (5.10). As discussed in Sec- 
tion 3, in the case of quantum cosmology one is not obviously obliged to take the same 
conditions, and indeed one can argue that, beyond (5.10), the boundary conditions are up 
for grabs. Indeed, one could effectively choose boundary conditions by proposing that the 
solution is given by the formula (5.1) with the factoring condition (5.10). As one can see 
from (5.12), this proposal again produces solutions with a time symmetric flavour to them 
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(i.e., with an equal number of G's and Gt's in the solution), a feature which persists to 
the case of more than one detector. Not surprisingly, this choice does in fact produce, in 
essence, the Hartle amplitude (3.8). 

We have, therefore, in this Section produced another candidate expression for the 
detection amplitude, which is arguably the more appropriate one for quantum cosmology. 
Furthermore, the means of generating it, Eqs.(5.1), (5.2), are readily generalizable to more 
complicated situations. 



6. AN IMPROVED DETECTOR MODEL 

We now briefly consider a more elaborate detector model that consists of a harmonic 
oscillator instead of the simple two state system. So we take, in the case of a single detector, 

Hfi = oja'^a, Hint = /(x)(a + J) (6.1) 

The energy eigenstate of the total Hamiltonian is again calculated using (5.2), with the 
factored fiducial state (5.10), where |0) is the harmonic oscillator ground state. In terms 
of a path integral, denoting the harmonic oscillator coordinates by q, 

*(x;,9y) = r dr e^^^ f V^{t)Vq{t) 

J —oo J 

X exp(i5o[x(t)]+i5^[g(t)]+iA^i„i[Q(t),x(t)]) Xo(xo)wo(?o) (6-2) 
where UQ^q) = {q\0). The integral over q is conveniently rewritten as 

exp {iS(i[q{t)] +iXSint[q{t),:>a{t)]) wo(go) 

= {qflTexp {-iT{Hd + XH^^t)) |0) (6.3) 
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where T denotes time ordering. The right-hand side of (6.3) is just the unitary evolution of 
the vacuum state for a driven harmonic oscillator. Using the properties of coherent states 
[38], (6.3) is equal to {qj\z{T)), where \z) is a standard coherent state and 

z{t) = -iX r dt e*'^(^-*) /(x(t)) (6.4) 

Jo 

To find the amplitudes for detection and no detection we expand the total state (6.2) in 
terms of the eigenstates of the harmonic oscillator. Since all states other than the ground 
state correspond to detection, there is no single amplitude corresponding to detection 
(although there is a probability). It is therefore easier to look at the amplitude for no 
detection, which is obtained by overlapping (6.2) with the ground state UQ{qf), yielding, 

/CO r 
dr e^^^ / Px(t) (0|^(r)) exp (i5o[x(t)]) xo(xo) (6.5) 
-oo J 

Prom the properties of coherent states [38] , we have 

(0|z(r))=exp(^-l|;^(r)|2) (6.6) 
The probability for no detection is {i^ndl'^nd) probability for detection is simply 

The result (6.5) clearly has the desired properties. For paths x(t) which never enter the 
detection region, z{t) = and the path integral is unaffected. Paths that enter the region, 
on the other hand, generally have z{t) ^ 0, and they are exponentially suppressed. This 
is therefore a much improved detector model in comparison to Eq.(5.8). Its validity is not 
restricted to the perturbative regime. The generalization to many detectors is trivial and 
essentially the same result concerning peaking about classical trajectories is the obtained. 

Of course, this detector is still not fully satisfactory because it can happen that z{t) = 
even for paths that enter the region, because of the oscillatory nature of z{t), hence we 
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again encounter the issue of detector recurrences. Again this will be avoided if the time 
the trajectory spends in the detector region is short (less than oj~^). 

A more challenging detector improvement avoiding the recurrence problem would be 
to construct one that, were it used in standard unitary quantum mechanics, would be 
irreversible, i.e., involves an essentially infinite number of degrees of freedom. Such a 
detector was introduced in the related context of measuring arrival times in Ref.[42] and 
it would be interesting to incorporate it into the situation considered here. 



7. SUMMARY 



The aim of this paper was to give substance to the appealing intuitive notion that so- 
lutions to the Wheeler-DeWitt equation (1.1) correspond to entire histories of the universe 
with time emerging as a parameter along each trajectory. The concrete technical results 
- the detection amplitude and the introduction of a set of timeless coherent states - are 
compatible with this notion. There are, however, many subtle aspects to this notion [4], 
and we do not claim to have an exhaustive demonstration of the emergence of trajectories 
from the Wheeler-DeWitt equation. 
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